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I. INTRODUCTION 

Anomalous transport phenomena occur in many phys- 
ical systems. Particles in complex media, instead of 
the normal diffusion they undergo anomalous diffusion. 
The basic property commonly used to characterize differ- 
ent types of diffusion is the mean-squared displacement 
(MSD) of the diffusing particles. The diffusion is known 
as anomalous diffusion if the MSD is no longer linear 
in time as in normal diffusion or Brownian motion; in- 
stead it satisfies a power law behavior and varies as t"' , 
< a < 2. MSD with a < 1 corresponds to subdiffusion, 
which is slower than the normal diffusion. It is known as 
superdiffusion when a > 1 with the diffusion faster than 
the normal one. Just like normal diffusion which can 
be described by diffusion equation and Langevin equa- 
tion, it is possible to model anomalous diffusion using 
fractional version of these equations. The stochastic pro- 
cesses commonly associated with anomalous diffusion in- 
clude fractional Brownian motion and Levy motion . 
In order to describe anomalous diffusion, the usual ran- 
dom walk model for normal diffusion has to be replaced 
by the continuous-time random walk model P, |3j HI • 

Anomalous diffusion forms only a portion of diffusion 
processes in nature. There also exist many diffusion pro- 
cesses which do not have MSD that varies as power law 
with a unique diffusion or scaling exponent a. For such 
systems a is no longer a constant, instead it can be a func- 
tion of certain physical parameters like position, time, 
temperature, density, etc. Kinetic equations of constant 
fractional order such as fractional diffusion equation and 
fractional Langevin equation are successful in describing 
anomalous diffusion. However, for diffusion process with 
non- unique scaling exponent, constant order fractional 
kinetic equations are not applicable. It is necessary to 
introduce kinetic equations of multi-fractional order or 



variable fractional order [6|-|8j We remark that variable 
order fractional differential equations are mathematically 
intractable except for few very simple cases. On the other 
hand, there is a subclass of diffusion processes which have 
different regimes with distinct scaling exponents. For ex- 
ample, single- file diffusion of Brownian particles confined 
to narrow channels or pores 0, |^, [13] • There also exists 
ultraslow diffusion such as the diffusion observed in Sinai 
model [ll| which has MSD that varies asymptotically as 
power of logarithm of time. It is possible to model such 
diffusion processes with a different type of multifractional 
differential equation based on distributed order derivative 
first introduced by Caputo [T2j . 

Since its introduction in the 1960s by Caputo [l^l, dis- 
tributed order differential e quat ion was subsequently de- 
veloped by Caputo himself |l3l . IT^ , and also by Bagley 
and Torvik [III, [l^. Various authors [iTl - fssj have ap- 
plied fractional differential equations of distributed or- 
der to model anomalous, non-exponential relaxation pro- 
cesses and anomalous diffusion with non-unique diffusion 
or scaling exponent. 

The distributed order time derivative is defined by 



D(^)f{t)= ^\ (a) D'^ fit) da 



(1) 



where the weight function ip{a) is a positive integrable 
function defined on [/3i,/32]. For our purpose, it is as- 
sumed that 0</3i<a</32<l- Here the fractional 
time derivative D" can be either the Riemann-Liouville 
or Caputo type [l], [s^, , which are defined respectively 
for n— l<a<nas 



D%Lfit) 



1 



f{u)du 



(2) 



and 
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T{n ~ a) Jq {t - u)" 



/(")(M)du 



(3) 



For our purpose, we let /3i = and /32 = 1 in the subse- 
quent discussion. 
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Due to the fact that distributed order derivative mod- 
ifies the fractional order of derivative by integrating all 
possible orders over a certain range, the solution of the 
resulting fractional equation is not characterized by a def- 
inite scaling exponent. One can regard the distributed or- 
der fractional time derivatives as time derivatives on var- 
ious time scales. Derivatives with distributed fractional 
order can be used to describe transport phenomena in 
complex heterogeneous media with multifractal property. 
Such processes exhibit memory effects over various time 
scales. Distributed order time-fractional diffusion equa- 
tion and space-fractional diffusion equation have been 
considered by various authors 0-[iS El-Ell- Solutions 
of these equations can be used to describe retarding sub- 
diffusion and accelerating subdiffusion, as well as acceler- 
ating supcrdiffusion fisl Eol . [22l . [2^ . [26l [3l| and ultraslow 
diffusion [HM ii, IH IM IM, l31| . 

Distributed order fractional diffusion equations of var- 
ious types have been quite well-studied. On the other 
hand, distributed order Langevin equations have not 
been considered as far as we are aware. The main aim of 
this paper is to study the Langevin-like equations of dis- 
tributed order and to consider their possible applications. 
This paper considers several types of fractional Langevin 
equation of distributed order. The statistical properties, 
in particular the MSD, of the solutions to these equations 
are studied. Possible applications of these equations to 
model retarding subdiffusion such as single-file diffusion 
and ultraslow diffusion such as in the Sinai model and 
some other systems will be discussed. 



II. FRACTIONAL LANGEVIN OF 
DISTRIBUTED ORDER 

Let us consider a simple free fractional Langevin equa- 
tion without frictional term 



and 



(4) 



where ^(t) is a stationary Gaussian random noise with 
mean zero and covariance C(t) = {x{t + T)x{t)) to 
be specified later on. Using the Laplace transform 
of Riemann-Liouville derivative D'^]^f(t) and Caputo 
derivative D^f{t) 0,[2l|: 



k=0 



^ [D%Lf{t)] {s) = - E [^Sl'"VW],=o (5) 

and 

C [D^fit)] (s) = sy~{s) - (0), (6) 



n-1 



fc=0 



one gets the Laplace transform of ([4]) for the Riemann- 
Liouville and Caputo cases are respectively 

n-1 
A:=0 



n-1 

E 

fc=0 



"i(s) - V s'^-^-^x'^'"'^ (0) = 



(8) 



The solution is 



1 



x{t)=ait) + -^ I {t-ur-'au)du (9) 

r(a) Jo 

where a(t) is the inverse Laplace transform of 

n— 1 n— 1 

^s''[D'^l''-'^x{t)]^^^ and ^s"-''^-ia;('=)(0) for the 

k=Q fe=0 

Riemann-Liouville and Caputo case respectively. If 
the random noise is white noise r](t) with covariance 
{'q{t)ri{s)) = d{t — s), then the variance of the process 
is given by 



1 



r(a)2 Jo 



t 



2a-l 



{2a - l)r(a)2 



, (10) 



and it is the same for both Caputo and Riemann-Liouville 
derivative. If one assumes a{t) = 0, ([9]) can be regarded 
as the definition of Riemann-Liouville fractional Brow- 
nian motion with the Hurst index H = a — ^ {sgI Wi\ . 
When a = 1, we get Brownian motion with H = 1/2. 

One can now generalize ^ to the case of distributed 
order fractional Langevin equation 



D^^^x{t)=^{t), t>0, 



(11) 



with -Df^) as defined as in ([T]). Note that though ([TT|) is 
in the form of free fractional Langevin equation of dis- 
tributed order, we shall see in next section that the dis- 
tributed order derivative term will introduce "frictional" 
terms depending on the type of weight function ^p{a). 
Laplace transform of (|lip is 



Ais)xis) - B{s) = as) 



(12) 



with 



A{s) = I ip{a)s°'da, < a < 1, (13) 
Jo 

for Riemann-Liouville case, and 

B{s) = da^{a) s^-'^-ia;^^) (0) ^ (14b) 

for Caputo case. Here \a] denotes the largest integer 
smaller or equal to a. 
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Solving (|T2)) gives 



i(s) 



B{s) i{s) 
A{s)^ A{sy 

with inverse Laplace transform 

x{t) = m{t) + [ G{t- T)^(r)dT, 



(15) 



(16) 



where m(t) and G{t) are the inverse Laplace transform 
of B{s)/A{s) and 1/^4(5) respectively. 

The mean and variance of the process x{t) are given 

by 



X = {x{t)) = rn{t), 



(17) 



and 



2 / G(u) / Ciu-v)Giv)dvdu 
Jo Jo 



(18) 



Thus, the solution to the distributed order fractional 
Langevin equation pi|) have the same variance for both 
the Riemann-Liouville and Caputo cases. Their means 
are different except when ^(''^(O) and [_D^^''"^a;(t)] = 0, 
fc = 0, 1, • • • ,n, and in this case the MSD for the two type 
of fractional derivatives are the same and equal to the 
variance. This is in contrast to time-fractional diffusion 
equation of distributed order, which leads to stochastic 
process with different MSD or variance for the Riemann- 
Liouville and Caputo case In our 
subsequent discussion D" and ^(,^) shall denote respec- 
tively fractional derivative and distributed order frac- 
tional derivative of either Riemann-Liouville or Caputo 
type. The two terms MSD and variance shall be used 
interchangeably 

Before we discuss various examples of distributed order 
Langevin equation, we first derive expressions of MSD 
for two types of Gaussian noise. First, we let ^(t) be the 
simple case of Gaussian white noise r]{t) with zero mean 
and covariance Cjj{t — s) = {i]{t)j]{s)) = 6{t — s). Then 
(ITSll becomes 



a^t) 



{G{u)ydu. 



(19) 



Evaluation of (jT9|) is in general complicated and it usually 
does not lead to a closed expression. One way to obtain 
a simpler expression for MSD is to impose the follow- 
ing condition on the Laplace transform of the covariance 
Q(s) of^(i): 



G{s)G{s) = -. 
The condition (gOl) reduces ^ to 

(j^{t) = 2 [ G{u)dv 
Jo 



(20) 



(21) 



As we shall show in subsequent sections that Gaussian 
random noise with Laplace transform of its covariance 
satisfying condition (|20p not only facilitate the calcula- 
tion of MSD, it also allows pTj) to model ultraslow diffu- 
sion which otherwise can not be done with the Gaussian 
white noise. 

Let us consider the properties of the Gaussian random 
noise ^(i) that satisfies condition (P0|) . From (|^ one 
gets 



C'is) 



1 



sG(s) 

= / 'rf5(a)s 
Jo 



Ms) 

s 

a-l. 



'da. (22) 

By noting that the inverse Laplace transform of s"^^ 
t-"/r(l - a), one gets 

f-i 



IS 



Q(t) 



t- 



-ip{a)dc 



< Of < 1. 



lo r(i-a) 

Let us define the Gaussian random noise S,a{t) by 



(23) 



and 

{Ut)Us)) 

such that 

cs-s)^m)m = 



{t-s) 



r(i 



(24) 



(25) 



Jo 



{ia{t)^p{s)) dadp 



{S.a{t)ia{s)) (p{a)da, 



(26) 



where we have used the orthogonality property of 
(Ca(i)C/3(s)) = Ha - /3) (^a(i)Ca(s))- Recall that the 
increment process of fractional Brownian motion or frac- 
tional Gaussian noise has the covariance Cnt^^^^ , where 
< -ff < 1 is the Hurst index for the fractional Brownian 
motion, and Gh = 2H{2H - 1). By letting a = 2-2H, 
one can then identify ^a(i) with the fractional Gaussian 
noise (up to a multiplicative constant). Note that the 
fractional Gaussian noise is to be regarded as generalized 
derivatives of fractional Brownian motion or a general- 
ized Gaussian random process [11] . We shall denote the 
distributed order fractional Gaussian noise with covari- 
ance (|26|) by ^-n^d it can be defined either as 



ip{a)^a{t)da, 



or 



(27a) 



(27b) 



depending on the nature of the weight function ip{a) as 
we shall demonstrate in subsequent sections. 
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Before we end this section, it is necessary to point out 
some problems associated with the stochastic differen- 
tial equation driven by fractional Gaussian noise. The 
question on whether a stochastic integral with respect 
to fractional Brownian motion leads to a well defined 
stochastic integral is a long standing problem which has 
attracted considerable attention (see [s^ HO] and refer- 
ences therein) . Fractional Brownian motion is not a semi- 
martingale if the Hurst index of the process H 7^ 1/2, 
that is when the process is not a Brownian motion. As a 
result, the usual stochastic calculus of Ito can not be 
used to define the integrals with respect to fractional 
Brownian motion. Various methods such as Sokorohod- 
Stratonovich stochastic integrals, Malliavin calculus, and 
pathwise stochastic calculus have been siiggcsted to over- 
come this problem (see references [sol . lift] for details). 
However, theory based on abstract integrals will en- 
counter difficulty in physical interpretations in certain 
applications. Since our subsequent discussion deals with 
applications involving 1/2 < H < 1, it is possible to 
consider the integrals with respect to fractional Brown- 
ian motion as the pathwise Ricmann-Stieltjes integrals 
(see for example j4l| and references given there). In this 
way we can handle such integrals in a similar manner as 
ordinary integrals. 



III. DOUBLE DELTA FUNCTION 
DISTRIBUTED ORDER FRACTIONAL 
LANGEVIN EQUATION 

Fractional Langevin equation of double order results 
with the choice of the following weight function: 

Lp{a) ~ ai^ia — oci) + a2(5(a — Q!2), (28) 

where 1/2 < a\ < a2 < 1. The distributed order 
Langevin equation ([TT]) becomes 

a2D°'^-x{t) +aiD'^'x{t) = ^{t) (29) 

For both Ricmann-Liouvillc and Caputo cases, one has 

A(s) = ais"i +a2s"^ (30) 

such that the Green function is given by the inverse 
Laplace transform of 



1 



1 



1 



A{s) ais"i-fa2s"2 02 s"2-"i -|- (01/02) ' 
That is. 



G(t^ — - t Ea2—Q-i,ct2 



0,2 



0-2 



(31) 



(32) 



where 



- r(/ij + i.) 



, n>0, v>0, (33) 



is the Mittag-LcfSer function [34 1. 

We consider first the case where the random noise is 
given by white noise ri{t), then the MSD of the process 
is given by 

a^{t)= f f G{u)5{u - v)G{v)dudv = f [G{u)f du 
Jo Jo Jo 



1 



'2 Jo 



u 



2(a2-l) 



aiu 



02 



du 
(34) 



which can not be evaluated analytically. However, its 
asymptotic limits can be obtained and are given by 



t 



2C(2 



a2(2a2-l)(r(a2))' 



and 



a\t) 



af(2ai-l)(r(ai))^ 



as t -> 



as i — cx). 



(35) 



(36) 



One can view the short time limit of MSD is obtained 
by ignoring the D°'^x{t) in ([29)l when i — > and treated 
like a2D°''^x{t) = ^(i). On the other hand, the long time 
limit of MSD is resulted when D°'^x{t) is neglected and 
p9)) reduces to aiD°'^x{t) = ^{t) as t ^ 00. Therefore 
the process described by (^5]) initially diffuses with scal- 
ing exponent a2 , then it slows down to become a process 
with a smaller scaling exponent ai as t becomes very 
large. In other words, the process is asymptotically lo- 
cally self-similar of order 012 — 1/2 with x{ct) = c^^'^xlt) 
for c > 0; and it is large-time asymptotically self-similar 
of order ai — 1/2 with x{ct) — c°''^~^x{t). Thus the re- 
sulting process describes retarding subdiffusion, which 
becomes more and more anomalous (or getting slower 
and slower) as time progresses. If we allow the lim- 
its of integration in (IT|) to be /3i = 1, /32 = 3/2, and 
1 < ai < a2 < 3/2 in ([28|) . then the resulting process 
is retarding superdiffusion. However, there is no way for 
one to obtain accelerating subdiffusion or superdiffusion 
based on (gH). 

Next we consider the case where the random noise ^{t) 
in pT|) is given by the distributed order fractional Gaus- 



sian noise 



S,ait)ip{a)da = V«iCai(i) + V^Ca2(0 (37) 



with covariance given by 



Q(i) =ai 



r 



0.2-, 



r(i-ai) -r(i-a2) 



(38) 



The random noise ^(t) can be regarded as the sum of two 
fractional Gaussian noises correspond to Hurst indices 
if, = l-f, 1 = 1,2 

Stochastic process associated with the distributed or- 
der fractional Langevin equation with the above 
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Gaussian random noise £,y^{t) has the following MSD: 



"2 Jo 



Ct2 —Cti ,Q2 



0-1 

1 

0-2 



du 



02 



Ol 

! 

02 



(39) 



The long and short time limits are then given by 
2 



and 



a2T{a2 + 1) 



as t -> (40) 



air(ai + 1) 



&Bt^oo (41) 



Thus we sec that the solution of (p9)) with Gaussian white 
noise and Gaussian distributed fractional noise both lead 
to power law type of MSD with different exponents for 
the short and long time limits of MSD. 

Now we want to see whether the processes in the above 
examples can be used to describe certain transport phe- 
nomena in physical systems. Clearly, such a process must 
not have unique characteristic scaling exponent, instead 
it has piecewise scaling exponents. One such physical 
process which has different scaling for short and long time 
behavior of MSD is single-file diffusion. Recall that in 
single-file diffusion, the particles are geometrically con- 
straint to move in a line and unable to overtake 0, 0, [13] ■ 
For very short times the MSD of the diffusing particles 
varies with time, just like in normal diffusion. The Brow- 
nian particles in confined geometry such as nanopores or 
nanochannels can not alter their relative ordering, there- 
fore the subsequent motion of each particle is always con- 
strained by the same two neighboring particles. Thus, in 
the long time limit this effect of caging slows down the 
diffusion, and changes the MSD from linear growth to 
one that varies with y/t. If white noise is used in p9|) . 
when a2 = 1 and 2ai — 1 = 1/2, or ai = 3/4, the limits 
(|35|) and ((36|) give the correct asymptotic behavior of the 
MSD for single-file diffusion [2,la[l3l- On the other hand, 
if we use distributed order Gaussian fractional noise in 
(P^ . then for a2 = 1 and ai = 1/2, one obtains the 
correct short and long time limits for the MSD of single- 
file diffusion. Figure 1(a) shows the comparison of MSD 
obtained using Gaussian white noise and distributed or- 



der fractional Gaussian noise, and Figures 1(b) and 1(c) 



show the short and long time limit of MSD for these two 
cases. 

Recall that Brownian motion is related to white noise 
(in the sense of generalized function) by the following free 
Langevin equation Dx{t) = ri{t). If we regard the pro- 
cess as Brownian motion experiences a retardation due 
to the confined geometry, then with Q!2 = 1 and ai = 3/4 
such that the second term with the fractional derivative 
in (|29p can be regarded as a damping or retarding term 
that slows down the Brownian motion so that the mo- 
tion which begins as normal diffusion becomes single-file 



(a) 




(b) 




FIG. 1. MSD of double delta-function distributed order pro- 
cess with Ol = 02 = 1, a2 = 1 and a\ — 3/4. (a) MSD associ- 
ated with distributed order equation with white noise (WN) 
and double delta function distributed order fractional Gaus- 
sian noise (FGNi); (b) log- log plot of MSD associated with 
double-order distributed order fractional Langevin equation 
with WN; (c) log-log plot of MSD associated with double- 
order distributed order fractional Langevin equation with 
FGNi. 



subdiffusion after a long time. Since we only use ([^^ to 
describe the asymptotic properties of single-file diffusion, 
it does not give a unique description of single-file diffu- 
sion. It is necessary to consider the behavior of single- 
file diffusion at intermediate times to see whether it also 
agrees with the corresponding description given by the 
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distributed order fractional Langevin equation (|29p . 

Here we would like to remark that a similar asymp- 
totic behavior for the MSD can also be obtained by us- 
in g fra ctional time diffusion equation of distributed order 



ip{a)D^W{x,t)da 



(42) 



where x e M, t > and W{x^ t) is the probability distri- 
bution function. Using the same weighing function ^p{a) 
as in and Z?" as Caputo fractional derivative one 
gets 



d'^^W{x,t) d°'^W{x,t) _ d'^W{x,t) 



02- 



dx'^ 



(43) 



Using the initial condition W{x^ 0+) ~ 5{x)^ the solution 
of (|43p is a diffusion process with variance having asymp- 
totic behavior tr^ ~ 1°^^ and ~ t"^ as large time and 
short time limit respectively. However, if the fractional 
time derivative of Riemann-Liouville type is used in (|42|l . 
then the asymptotic behavior is opposite to that for the 
Caputo derivative with the smaller exponent ai domi- 
nates at small time and larger exponent a2 dominates at 
large time. Thus, in contrast to ((29|) . it is possible for to 
obtain accelerating subdiffusion based Riemann-Liouville 
version of (|42p . Another major difference between these 
two approaches is that the description based on the dis- 
tributed order time fractional diffusion equation (|42|) is a 
non-Gaussian model, whereas the distributed order frac- 
tional Langevin equation pQ]) is a Gaussian one. We 
would like to remark that there also exists an effective 
Fokker-Planck equation which leads to a similar results, 
and it provides a Gaussian model for single-file diffusion 

El. 



IV. UNIFORMLY DISTRIBUTED ORDER 
FRACTIONAL LANGEVIN EQUATION 

There exists a class of strongly anomalous diffusion 
with the long time limit of its MSD decays logarithmi- 
cally as {\ogtY\ K > 0. Such ultraslow diffusion oc- 
curs in Sinai diffusion of a particle in a one-dimensional 
quenched random energy landscape [ll|, , in charged 
polymers (43 |. motion in aperiodic environments [45j . 
in a class of iterated maps [46j . in area preserving 
parabolic map [l^ and charged tracer particle on a two- 
dimension lattice [i^, etc. It has been shown that uni- 
formly distributed order time fractional diffusion equa- 
tion can be used to model the ultraslow diffusion [isl - 
[2ll,[2i,[2l,[2l[3S[3l|. 

In this section we want to consider uniformly dis- 
tributed order fractional Langevin equation to see 
whether it can describe ultraslow diffusion. For uni- 
form distributed order the weight function is ^p{a) — I, 



< a < 1. Now (fTTj) becomes 



D°'x{t)da = ^{t) 



(44) 







which gives the Laplace transform of its Green function 
as 



G(,s) 



1 



A{s) 



s"da 



s - I 



(45) 



By taking the inverse Laplace transform one gets 

G(t) =e*Ei(i), (46) 

with 

Ei(t) = -7-logi + Ein(t), (47) 
Ei(t) is the exponential integral function given by [i^ 



Ei(z) 



y 



-dy 



(48) 



and 



Ein(t) 



I - e-'" 



E^^rr^- (49) 



fe=i 



kkl 



In the case where the random noise ^ (t) is white noise, 
the MSD is given by 



a^{t) = f [e"Ei(u)]^dM, 
Jo 



(50) 



which can not be evaluated analytically. In order to study 
the asymptotic behavior of the MSD, we consider the up- 
per and lower bound of the exponential integral function 
(see m, #5.1.20): 



^\og(l + j\ <e*Ei(0 <log + ^ 



From (|50p the upper bound of the variance is 



(51) 



U{t) 



log 1 + - 

u 



3 



+ t 



loa 



1 + t 



du 

2 



oo ^ 



^^n2(i + t) 



(52) 



The lower bound is 





2 70 



log 1 + - 
u 



log 1 



du =-u(- 
2 I 2 



(53) 



When t — > 0, the summation term in ([52|l tends to zeta 
function 2^(2) = tt^/B, which cancels with the first term 
in the equation. Therefore one gets 



U{t) t [log(t)] 



(54) 
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and 



Thus, the short-time hniit of the MSD is given by 



(55) 



<J^{t) cit [log(i)]' = cit [log(lA)]' (56) 

where 1/4 < ci < 1. From numerical simulations, we 
get ci = 1 as shown in Figure |2(b)[ For the lar ge-time 
limit of the MSD, one notes that the summation term in 
([52]) tends to zero as t — > oo, and the second term of ((52|) 
becomes 







1 


log 


= t 









1 

2fi 



Thus, 



and 



U{t) 



TT 

t—^oo 3 



(57) 



(58) 



(59) 



Therefore the MSD approaches a constant for sufficient 
large time. 



o'{t) ^ C2 — 



(60) 



where C2 = 1.5 can be obtained graphically. One may 
retain the time-dependent term in the MSD for large (but 
not too large) time: 



2. ^ 1 

{t) - -r + 7 

t^oo 4 t 



(61) 



Here we have a motion that begins as a non-stationary 
process and becomes a stationary one after sufficiently 
long time. In fact, at small times the diffusion is anoma- 
lous diffusion of slightly superdiffusive type. In other 
words, the process goes to zero slower than the normal 
diffusion due to the log(l/i) term as < — >■ 0. However, at 
large times it tends to a stationary process with a con- 
stant variance. One can interprete the long time behavior 
in the following way. The uniformly distributed deriva- 
tive is the derivative D" integrated over the range a = 
to a = 1. As i — >■ cxj, one would expect the dominant 
term will be from a = 0, which will result x{t) ~ illt), a 
white noise process with constant variance. In compari- 
son, for distributed order time fractional diffusion equa- 
tion (|^^ with uniformly distributed order, the associated 
process at short time behaves somewhat superdiffusive, 
with variance ~ tlog{l/t); and for long time limit the dif- 
fusion process becomes ultraslow with variance ~ log(t). 



Now instead of using white noise in (|44p . we let the 
random noise ^(t) be the uniformly distributed fractional 
Gaussian noise, 



S,ait)da, 



(62) 



where the weight function is given by ip{a) — 1. From 
([5^ one gets the covariance of (t) as 



C^{t - s) = {Ut)U^)) ^ 

The MSD is then given by 



r(i-a) 



da. 



(63) 



a^it)^2 G{u)du^2 e"Ei(u)du 
Jq Jq 

-2(e*Ei(t)+7 + log(t)) 



(64) 



where 7 is the Euler's constant. Substituting P7)) into 
(|64)) gives 

<7^ {t) = 2e*Ein(t) + 2 (l - e*) (7 + log(t)) (65) 
Hence, the short time limit of the MSD is given by 

a^(t)^^^2tlog(i). (66) 

Using ()5ip in ()64p . one gets the upper bound and lower 
bounds of MSD as 



H{t) = 2 



loa 



l + t 
t 



+ 7 + log(0 



= 2[log(l+i)+7] 
2 + t 



- loa 
2 



t 



- 7 + log(t) 



log [t(2 + t)] + 7 



(67a) 



(67b) 



Both upper and lower bound are approach to the same 
assymptotic function, thus we have 

a'^it) - 2\og{t)+2j (68) 

t-^OO 

Therefore, the large time hmit of the MSD is 

<j^(t) ^ 21og(t) (69) 

Note the constant term 7 « 0.57722 can not be neglected 
in practice since even for t = e^", a'^{t) « 2(10-1-0.57722) 
which shows a 6 % contribution from the Euler's con- 
stant. The above asymptotic behavior of the MSD shows 
that the diffusion is ultraslow at very large times, and it 
becomes slightly superdiffusive at small times. 

Figure 2(a) shows comparison of the MSD of the 



stochastic process associated with uniformly distributed 
order Langevin equations with Gaussian white noise and 
uniformly distributed fractional Gaussian noise. In Fig- 
ures [2(b)] and 2(c) the short and long time limit of MSD 
for these two cases is demonstrated. 
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(a) 




(b) 




FIG. 2. MSD of the uniformly distributed process corresponds 
to the uniformly distributed order in the range a £ [0, 1]: 
(a) MSD of the process associated with distributed uni- 
formly order fractional Langevin equation with white noise 
(WN) and uniformly distributed fractional Gaussian noise 
(FGN2); (b) log- log plot of MSD associated with uniformly 
distributed order fractional Langevin equation with WN; (c) 
log- log plot of MSD associated with the distributed order frac- 
tional Langevin equation with FGN2. 



V. POWER LAW DISTRIBUTED ORDER 
FRACTIONAL LANGEVIN EQUATION 

In order to describe ultraslow diffusion processes with 
large time limit of MSD varies as (logt)'^, > 0, it is nec- 
essary to consider power-law distributed order fractional 



Langevin equation with weighing function vo" ^, > 0. 
(fTTjl now becomes 



'0 



(70) 



Laplace transform of the Green function G{s) of ([70]) is 
the inverse of 



A[s) 



via 
la 



"da. 



For < s < 1 or —00 < log s < 0, one obtains 

(-logs)"^ 



1^7 (z^, -logs)' 



(71) 



(72) 



Since ^(v, z) ^ r(i/) as z — >■ 00, thus for small s or large 
- log s, 



G{s) 



(log(l/^))'^ 



as s — > 0. 



(73) 



If we assume the Gaussian random noise in (j70p is power- 
law distributed order Gaussian fractional noise, 



Jo 



(74) 



then the Laplace transform of the MSD is given by 
CT^(s) = 2G{s)Js, which can be verified as a slowly vary- 
ing function [sO, • Now applying Tauberian theorem 
[51I . [52! , which allows the long and short time asymptotic 
limits of a function f{t) to be obtained from the Laplace 
transform /(s) for s near origin and infinity respectively 
[see for example, page 445, reference 5l| ). Thus from 
([73]) . one gets 



'Tiv + l) 



a.s t 00. 



(75) 



Similarly, the short time limit for the MSD can be ob- 
tained from the large s limit of ct^(s) given by 



logs 



as s — >■ 00. 



(76) 



From this we get the short time limit for the MSD as 



tlog{l/t) 



as t 



0. 



(77) 



Thus, the distributed order fractional Langevin equation 
with power law weight function provides a way to de- 
scribe the kinetics of the ultraslow diffusion such as Sinai 
diffusion with v ~ A for particles moves in a quenched 
random field, and transport of hooked polyampholytes 
(heteropolymers which carry both positive and negative 
charges) with ly = 4/3. 
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VI. CONCLUDING REMARKS 

We have shown that distributed order fractional 
Langevin equation provides a mathematical model for 
anomalous diffusion which does not have a unique scal- 
ing exponent. It is interesting to note that the expres- 
sion for MSD acquires a more simple form if white noise 
in the distributed order fractional Langevin equation is 
replaced by distributed order fractional Gaussian noise. 
The solutions of distributed order fractional Langevin 
equations have MSD which describe retarding subdifFu- 
sion such as in single-file diffusion, and ultraslow diffu- 
sion with logarithmic growth. To a large extent, the re- 
sults obtained are similar to that from time-fractional 
diffusion equation of distributed order, except for one 
main difference that the MSD for the Langevin case has 
same properties for both Riemann-Liouville and Caputo 
distributed derivatives, whereas in the fractional diffu- 
sion equation Riemann-Liouville and Caputo distributed 
derivatives lead to MSD with different behavior. In addi- 
tion, distributed order time fractional diffusion equation 
result in a non-Gaussian process, whereas in the process 
obtained from the corresponding Langevin equation is 
Gaussian. 

Possible direct generalizations of our study are exten- 
sions of free Langevin equation to fractional Langevin 
equation and fractional generalized Langevin equation 
of distributed order. However, one notes that frictional 
terms appear in the free Langevin equation once the usual 



fractional derivative is replaced by the distributed or- 
der fractional derivative. For example, free fractional 
Langevin equation of distributed order with the weight 
hmction ip{a) = aiS{a — ai) + a2(5(a — a2) + XS{a), 
A > OwiU become a2D°''x{t) + aiD"^x{t) + Xx{t) = ^{t). 
Hence frictional terms can appear in free Langevin equa- 
tion when the time derivative is replaced by distributed 
order time derivative with certain weight function. The 
solution is complex in such a case as the Green function 
involves sum of Wright functions [s^, and the associ- 
ated MSD is even be more complicated to compute. In 
the case of generalized fractional Langevin equation of 
distributed order D^^^xit) + j{t — u)x{u)du = ^(t), 
where j{t) is the frictional kernel, one would expect 
it to be mathematically more involved, though the as- 
sumption of fluctuation-dissipation theorem may help to 
simplify the situation somewhat. The main difficulty is 
in the evaluation of the inverse Laplace transform for 

the Green function G{s) = G^{s) 1 + G ,^,{3)^ [s) 

where G^{s) denotes Laplace transform of Green func- 
tion for D(^^p)x{t) = 0. All these generalizations are not 
only computationally complex and mathematically in- 
tractable, they may not lead to very interesting results. 
Our study shows that the simpler fractional Langevin 
equation of distributed order seem to be adequate for 
describing the kinetics of the types of diffusion under con- 
sideration, hence it provides a viable alternative to the 
time fractional diffusion equation of distributed order. 
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